Observation of metastable states in spinor Bose-Einstein condensates 
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Bose-Einstein condensates have been prepared in long-lived metastable excited states. Two comple- 
mentary types of metastable states were observed. The first is due to the immiscibility of multiple 
components in the condensate, and the second to local suppression of spin-relaxation collisions. Re- 
laxation via re-condensation of non-condensed atoms, spin relaxation, and quantum tunneling was 
observed. These experiments were done with F = 1 spinor Bose-Einstein condensates of sodium 
confined in an optical dipole trap. 
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Metastable states of matter, excited states which relax 
only slowly to the ground state, are commonly encoun- 
tered. This slow relaxation often arises from the pres- 
ence of free-energy barriers that prevent a system from 
directly evolving toward its ground state; if the ther- 
mal energy to overcome this barrier is not available, the 
metastable state may be long-lived. 

Many properties of Bose-Einstein condensates in di- 
lute atomic gases arise from metastability; indeed, 
such condensates are themselves metastable, since the 
true equilibrium state is a solid at these low temper- 
atures. Bose-Einstein condensates in gases with at- 
tractive interactions (scattering length a < 0) || are 
metastable against collapse due to a kinetic energy bar- 
rier H . The persistence of rotations in condensates with 
repulsive interactions (a > 0) hinges on whether vortices 
are metastable in singly- || or multiply-connected 
geometries. Similarly, dark solitons in restricted ge- 
ometries are predicted to be long-lived akin to the 
recently observed metastable states in superfluid 3 He- 
B pj. Finally, Pu and Bigelow discussed spatial distri- 
butions of two-species condensates which are metastable 
due to mean- field repulsion between the two species [[ll] . 

In this Letter, we report on the observation of two 
complementary types of metastability in F = 1 spinor 
Bose-Einstein condensates of sodium. In one, a two- 
component condensate in the \F — l,mp = 1,0) hyper- 
fine states was stable in spin composition, but sponta- 
neously formed a metastable spatial arrangement of spin 
domains. In the other, a single component \mp = 0) con- 
densate was metastable in spin composition with respect 
to the development of \mp = ±1) ground-state spin do- 
mains. In both cases, the energy barriers which caused 
the metastability (as low as 0.1 nK) were much smaller 
than the temperature of the gas (about 100 nK) which 
would suggest a rapid thermal relaxation. However, this 
thermal relaxation was slowed considerably due to the 
high condensate fraction and the extreme diluteness of 
the non-condensed cloud. 

Spinor Bosc-Einstcin condensates were created as in 



previous work |12|. Condensates in the \F = l,mp = 
— 1) state were created in a magnetic trap fD| and then 
transfered to an optical dipole trap formed by a single 
infrared laser beam [[u} . The beam was weakly focused, 
producing cigar-shaped traps with depths of 1 - 2 /xK, 
radial trap frequencies of about 500 Hz, and aspect ra- 
tios of about 70. After transfer, Landau-Zener rf-sweeps 
placed the optically trapped atoms into the desired initial 
spin states. Condensates were then held for as long as 
20 s during which time the condensate numbers decayed 
from about 10 6 to 10 5 . Condensate densities ranged from 
about 5 to 1 x 10 14 cm~ 3 . 

The spin-state distribution along the long axis of the 
optical trap was probed by a Stern-Gerlach spin separa- 
tion combined with time-of-flight imaging Q. After a 
sudden switch-off of the optical trap, the condensate ex- 
panded primarily radially. After 5 ms, a magnetic field 
gradient of several G/cm separated the different spin- 
state populations without distorting them. Axial stria- 
tions observed in time-of-flight images therefore reflected 
the spatial structure of the trapped condensate. After a 
total time-of-flight of 15 - 30 ms, the atoms were opti- 
cally pumped to the \F = 2,mp — 2) hyperfine state, 
and imaged with resonant light. 

For our first set of experiments, a condensate in a su- 
perposition of the \F = l,mp = 1,0) hyperfine states 
was prepared by first transferring all the atoms to the 
\mp = 1) state with an adiabatic rf-sweep at low field 
(1 G), and then transferring half of the population to 
the \mp = 0) state with a non-adiabatic sweep at high 
field (15 G). The sample then evolved freely in a 15 G 
bias field which prevented the formation of \mp = — 1) 
spin domains due to the large quadratic Zeeman energy 
barrier [fiJJ. 

Following their preparation, the two spin components 
began to separate radially (Fig. 1, sideways in images). 
Within 50 ms, axial density striations were observed, in- 
dicating the formation of \mp = 1) and \mp — 0) spin 
domains. The striations were initially angled due to ra- 
dial excitations which tilted the domain walls at various 
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angles to the long axis of the cloud. After about 100 ms, 
the striations became horizontal, indicating the dissipa- 
tion of radial motion and the orientation of the domain 
walls perpendicular to the axis of the cloud. The spin do- 
mains reached typical sizes of 40 ± 15 /im, slightly larger 
than their size after 50 ms. Thereafter, the clouds were 
essentially unchanged - only the total number of atoms 
diminished due to three-body losses. Matching the im- 
ages of the \mp = 1) and \itif — 0) spin distributions 
showed that these long-lived states consist of persistent, 
alternating \mp — 1) and \mp = 0) spin domains. 

Phase-separation in two-component condensate mix- 
tures is predicted by mean- field theory $J^L§. The 
mean-field interaction energy of such condensates is given 
by 2nh 2 /m x {n 2 a a a + n 2 ab + 2n a riba a b) where m is the 
common atomic mass, n a and rib are the densities of 
each of the components, a a and a& are the same-species 
scattering lengths, and a a b is the scattering length for 
interspecies collisions. The mean-field energy is mini- 
mized by phase-separation when the scattering lengths 
obey the relation a a b > y/a a ab- At low magnetic fields, 
the scattering lengths in the F — 1 spinor system 
are determined by ap tot= o and a,F tot =2 which describe 
collision between atoms with total angular momentum 
F tot = 0,2 §]]. Defining a = {2a Ftot=2 + a Ftot=0 )/3 
and Aa = (ai? tot =2 — ap tot= o)/3, the scattering lengths 
for the \mp — 1,0) two-component system are given 
by ao = a, and a\ = aoi = a + Aa. Since Aa was 
measured to be 3.5 ± 1.5bohr > 12|, aoi > ^Ja^ai 
and the components phase-separate. Interestingly, this 
phase-separation should not occur in the non-condensed 
cloud because the same-species mean-field interaction en- 
ergies are doubled due to exchange terms. 

The tendency of a two-component mixture to phase- 
separate implies the presence of imaginary frequencies 
for out-of-phase collective excitations of the overlapping 
components (spin waves) [^9|,|l|. This means that the 
condensates are unstable, and slight perturbations grow 
exponentially. In our experiment, density and magnetic 
field inhomogeneities may yield slight perturbations in 
the initial state, corresponding to an initial distribution 
of spin wave excitations. The spontaneous formation of 
spin domains thus constitute an observation of spin waves 
with imaginary frequencies at wavelengths of 10 - 50 /im, 
corresponding to the size of the spin domains which are 
intially formed. 

The many-domain spin distribution is a macroscopi- 
cally occupied excited state; the ground state of the two- 
component system would contain only one domain each 
of \rriF = 0) and \mp = 1) atoms, minimizing the surface 
energy of the domain walls |24| . For the many-domain 
state to decay directly to the ground state, the two con- 
densate components would have to either overlap, or else 
pass by each other without overlapping. Such motion 
is energetically forbidden: the former due to the mean- 
field energy (typically 50 Hz or 2.5 nK |l2|), and the 



latter due to the large kinetic energy required to vary 
the condensate wavefunction radially (the condensate is 
about 5 /im wide, yielding a kinetic energy barrier of 
50 - 100 Hz). For this reason, the excited states we 
observe are metastable. In comparison, experiments on 
two-component mixtures of 87 Rb showed no metastabil- 
ity, due in part to the nearly spherical traps (aspect ratio 
of 1/V8) which were used [E2|. 

The typical size of the metastable spin domains can be 
understood by considering the width of the domain walls 
(~ 3 /im [p3|). Spin domains smaller than this width do 
not prevent the tunneling of unlike-spin atoms through 
the domain, and thus small spin domains are free to mi- 
grate through the cloud and coalesce with other like-spin 
domains. Spin domains which are several times larger 
than the boundary width prevent the tunneling of atoms 
across the domain, and are thus stable. Metastable do- 
mains shorter than w 25 /im were not observed, although 
the above description suggests that domains as short as 
w 10 /im could be stable. This may be due to either 
an axial smearing of small domains during their free ex- 
pansion prior to imaging, or perhaps the destabilization 
of small domains by the initial radial excitations in the 
trapped cloud. 

The spontaneously formed metastable states were ex- 
tremely long lived. In the absence of an externally ap- 
plied magnetic field gradient, the states persisted for at 
least 20 seconds, by which time the condensate num- 
ber had dropped so low that we were unable to prop- 
erly assess the time-of-flight images. The lifetime of the 
metastable state was shortened by applying a magnetic 
field gradient along the axis of the trap. At field gradients 
of 0.1 G/cm or greater, the metastable states decayed to 
the ground state within experimentally accessible times 
of 10 seconds or less. 

To study this decay with a well-defined initial state, we 
first created a metastable two-domain system. A strong 
magnetic field gradient was applied, causing the many 
spin domains to collapse into the ground state of two do- 
mains with the Imp = 1) domain at the high-field end 
of the cloud. A weaker field gradient (0.1 G/cm) was 
then applied in the opposite direction, which energeti- 
cally favored the rearragement of the spin domains at 
opposite ends of the cloud. In spite of this gradient, the 
two-domain system was metastable (Fig. 2). It decayed 
slowly to equilibrium via two processes: a slow process 
which dominated for the first 12 s, and a fast process 
(between 12 s and 13 s) which suddenly depleted the 
metastable domains. 

The slow process caused a gradual increase in the num- 
ber of atoms in the two ground state spin domains, which 
saturated at about 5 x 10 4 atoms in each. The rate of ac- 
cumulation did not depend strongly on condensate den- 
sity, which was varied by decreasing the total number 
of atoms before creating the metastable state. Thus, the 
slow process is probably due to re-condensation of the di- 
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lute thermal cloud into ground state spin domains, driven 
by a difference Afj, between the chemical potential of the 
thermal cloud and the energy of the ground state. 

Studies of the formation of one-component Bose con- 
densates in harmonic traps |^5],^6| considered the case 
T>p where T is the temperature of the cloud, and fi 
is the chemical potential of the condensate (fi > due 
to interactions). For optically trapped Bose condensates, 
the temperature is about 1/10 of the optical trap depth 
due to evaporation, i.e. T ps 100 nK puj ]. For typical 
densities of 3 x 10 14 cm -3 , fj, = 220 nK > T. Thus, pre- 
vious studies of condensate formation may not apply to 
these metastable states. 

The sudden decay to the ground state after the slow 
relaxation indicated the onset of a faster relaxation mech- 
anism, which occurred once the decaying condensate had 
reached a critical density. The dependence of this critical 
density on the applied field gradient suggests the onset 
of macroscopic quantum tunneling across the metastable 
spin domains. Such behaviour will be reported in greater 
detail elsewhere. 

The metastable states discussed so far involved 
multi-component mixtures stable in composition, but 
metastable in their spatial distribution. Conversely, we 
identified another form of metastability in which a sin- 
gle component condensate was stable in its spatial dis- 
tribution, but metastable in its hyperfine composition. 
This was observed in spinor condensates with an over- 
all spin projection (F z ) = along the magnetic field 
axis |l2] |. Condensates were prepared using Landau- 
Zener rf-sweeps at high fields (30 G) in one of two ways: 
either the entire trapped condensate was placed in the 
\mp = 0) state by a single rf-sweep, or the condensate 
was placed in a 50-50 mixture of the \mp = +1) and 
\m,F — — 1) states using two rf-sweeps. The condensates 
were then allowed to reach the ground state at a magnetic 
bias field Bq and axial field gradient B' . 

The evolution to equilibrium from the two starting con- 
ditions was quite different (Fig. 3). When starting from 
the \mp = 0) state, the condensate remained in that state 
for about three seconds before evolving over the next few 
seconds to the ground state. On the other hand, starting 
from an equal mixture of \mp = +1) and \mp = — 1) 
atoms, the fraction of atoms in the \mp = 0) state grew 
without delay, arriving at equilibrium within just 200 ms. 

This difference can be understood by considering a 
spin-relaxation collision, in which two |tof = 0) atoms 
collide to produce a \mp = 1) and a \rriF = —1) atom. 
In the presence of a magnetic field Bq, quadratic Zee- 
man shifts cause the energy of the two \mp = 0) atoms 
to be lower than the \mp = 1) and \mp = —1) atoms 
by 2 x 390 Hz x (B /G) 2 . Due to this activation energy, 
condensate atoms in the \uip = 0) state cannot undergo 
spin-relaxation collisions. Thus, even though the creation 
of \mp = 1) and \mp = —1) spin domains at the ends 
of the condensate is energetically favored globally in the 



presence of a magnetic field gradient, the \mp = 0) con- 
densate cannot overcome the local energy barrier for spin- 
relaxation. On the other hand, condensate atoms in the 
\mp = 1) and \mp = —1) states can directly lower their 
energy through such collisions, and equilibrate quickly. 

The decay of the metastable \mp = 0) state was stud- 
ied by varying the bias field Bq, the field gradient B 1 , 
and the trap depth U . Evolution curves under various 
conditions (Fig. 4) were characterized by a delay time Td 
before the fraction of condensate atoms in the \nip = 0) 
state decayed exponentially to equilibrium with a relax- 
ation time constant t. Increasing the bias field Bq did not 
affect r significantly, but increased the delay time from 
Td — 1.9 s at small bias fields (55 mG) to Td ~ 3.6 s for 
larger fields (250 mG). Lowering the trap depth from 2 
/iK to 1.2 /iK increased both Td from about 0.6 to 4 s, 
and r from about 3.8 to 6.6 s. Increasing the magnetic 
field gradient only induced a weak lowering of r, and Td 
was not affected significantly. 

The delay time Td should reflect the time-scale for the 
gradual accumulation of \mp — ±1) atoms in the thermal 
cloud before they reach a critical density and condense 
into the ground-state spin domains. Increasing the bias 
field increases the quadratic Zeeman energy barrier for 
spin-relaxation, lengthening Td- However, the clear vari- 
ation of Td as the energy barrier is changed from 0.1 nK at 
55 mG to 2.5 nK at 250 mG, energies much smaller than 
the temperature of T w 100 nK, cannot be explained on 
the basis of thermal spin- relaxation collisions. Further 
study of these relaxation timescales is warranted. 

In conclusion, we have identified a novel form of 
metastability in Bose-Einstein condensates which occurs 
at temperatures much higher than the energy barrier for 
evolution to the ground state. This metastability arises 
due to the high condensate fraction which dramatically 
slows down thermal relaxation. This offers the opportu- 
nity to study the kinetics of condensation and of the dis- 
solution of metastable condensates in "slow motion," on 
times scales much longer than collision times and of pe- 
riods of collective excitations. Furthermore, we have ob- 
served spontaneous phase-separation in a two-component 
mixture in accordance with the predictions of mean-field 
theory. 
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FIG. 1. Spontaneous formation of metastable states in a 15 
G bias field. Condensates were probed at various times after 
preparing an overlapping |?71f = 1,0) mixed condensate. The 
\rtiF = 1) (top) and \rtiF = 0) (bottom) states were separated 
during their free expansion. Striations reflect the presence 
of alternating ~ 40 um spin domains which form 50 ms after 
state preparation. High resolution (w 10 um) was achieved 
by only optically pumping a 200 um slice of the clouds. The 
height of each image is 1.3 mm. 
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FIG. 2. Relaxation of metastable states in a 15 G bias field 
and a 0.1 G/cm gradient. Images show metastable (outer) and 
ground state (inner) spin domains of \m,F = 1) and \m,F — 0) 
atoms, probed at various times after state preparation. The 
metastable domains decayed slowly for 12 s before tunneling 
quickly to the ground state. The height of each image is 1.3 
mm. 
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FIG. 3. Metastability of the pure \rriF = 0) state in the 
presence of a magnetic bias field (250 mG), and gradient 
(44 mG/cm). The evolution toward equilibrium of a pure 
\rriF ~ 0) condensate (open symbols), and a mixture of 
\m,F = 1) and \ttif = —1) (closed symbols) is shown by plot- 
ting the observed fraction of atoms in the \rriF — 0) state vs. 
dwell time in the optical trap. 
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FIG. 4. Decay of \rriF — 0) metastable states. The lines 
are exponential decay fits to the \mp = 0) fraction after a 
delay time Td- Increasing the bias field (top) increased Td- 
Lowering the temperature of the gas (bottom) increased both 
Td and r. The magnetic bias field Bo, the gradient B' , and 
the trap depth U for each curve are: a) 55 mG, 44 mG/cm, 
2 uK; b) 250 mG, 44 mG/cm, 2 uK; c) 20 mG, 11 mG/cm, 2 
uK; and d) 20 mG, 11 mG/cm, 1.2 uK. 
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